Multivariate characterisation of morpho-biometric traits of indigenous helmeted Guinea fowl (Numida meleagris) in Nigeria

This study was conducted to characterise phenotypically helmeted Guinea fowls in three agro-ecologies in Nigeria using multivariate approach. Eighteen biometric characters, four morphological indices and eleven qualitative physical traits were investigated in a total of 569 adult birds (158 males and 411 females). Descriptive statistics, non-parametric Kruskal–Wallis H test followed by the Mann–Whitney U and Dunn-Bonferroni tests for post hoc, Multiple Correspondence Analysis (MCA), Univariate Analysis, Canonical Discriminant Analysis, Categorical Principal Component Analysis and Decision Trees were employed to discern the effects of agro-ecological zone and sex on the morphostructural parameters. Agro-ecology had significant effect (P<0.05; P<0.01) on all the colour traits. In general, the most frequently observed colour phenotype of Guinea fowl had pearl plumage colour (54.0%), pale red skin colour (94.2%), black shank colour (68.7%), brown eye colour (49.7%), white earlobe colour (54.8%) and brown helmet colour (72.6%). The frequencies of helmet shape and wattle size were significantly influenced (P<0.01) by agro-ecology and sex. Overall, birds from the Southern Guinea Savanna zone had significantly higher values (P<0.05) for most biometric traits compared to their Sudano-Sahelian and Tropical Rainforest counterparts. They were also more compact (120.00 vs. 110.00 vs. 107.69) but had lesser condition index (7.66 vs. 9.45 vs. 9.30) and lower long-leggedness (19.71 vs. 19.23 vs. 9.51) than their counterparts from the two other zones. Sexual dimorphism (P<0.05) was in favour of male birds especially those in Southern Guinea Savanna and Sudano-Sahelian zones. However, the MCA and discriminant analysis revealed considerable intermingling of the qualitative physical traits, biometric traits and body indices especially between the Sudano-Sahelian and Tropical Rainforest birds. In spite of the high level of genetic admixture, the Guinea fowl populations could to a relative extent be distinguished using wing length, body length and eye colour. Generally, the birds from the three zones appeared to be more homogeneous than heterogeneous in nature. However, further complementary work on genomics will guide future selection and breeding programs geared towards improving the productivity, survival and environmental adaptation of indigenous helmeted Guinea fowls in the tropics.

breeding programs, to revaluate local breeds, allow the preservation of animal biodiversity and support consumer demands [22,23]. When such morphometric traits are considered jointly, multifactorial analyses have been shown to assess better the within-population variation which can be utilized in the discrimination of different population types [22,24].
In Nigeria, south Saharan Africa, there is dearth of information on the phenotypic diversity of Guinea fowls [25]. The current study aimed to find differences in indigenous Guinea fowl based on qualitative physical traits, biometric traits and morphological indices in three agroecological zones in Nigeria. The knowledge of the morpho-biometrical traits will support the implementation of breeding and conservation strategies in order to guarantee the survival and continuous production of the Guinea fowl genetic resource in the tropics for improved food security and livelihoods.

Ethics statement
In order to properly carry out the research, we adhered strictly to the ethical guidelines of the global code of conduct for research in resource-poor settings [26], following the Convention on Biological Diversity and Declaration of Helsinki. Although the study did not involve collection of blood and other tissue samples, we obtained field approval from the Research and Publication Directorate of Nasarawa State University, Keffi through permit no NSUK/FAC/ANS/ GF100. Written informed consent was also obtained from each participating farmer in line with best global practices.

Study area
The study was carried out in three agro-ecological zones of Nigeria namely; Sudano-Sahelian zone (Bauchi and Kano States), Southern Guinea Savanna zone (Nasarawa State and Abuja) and Tropical Rainforest zone (Ogun and Oyo States) (Fig 1). The Sudano-Sahelian zone is located between latitudes 10˚N and 14˚N and longitudes 4˚E and 14˚E, and lies immediately to the south of Sahara desert. The rainfall in this zone is less than 1000 mm per annum [27]. Temperature is high throughout the year with a mean minimum value of about 23˚C and mean maximum of about 34˚C. The zone is characterized by semi-arid grasslands vegetation while the density of trees and other plants decrease as one moves northwards. The Southern Guinea Savanna (GS) is part of the wider GS zone found on latitudes 7˚and 10˚N and longitudes 3˚and 14˚E [28]. The average annual maximum temperature ranges from 31 to 35˚C while the average annual minimum temperature is between 20 and 23˚C. It has mean annual rainfall of at least 1,600 millimeters and lowest mean monthly relative humidity of not less than 70 percent. It is a belt of mixture of trees and tall grasses. The Tropical Rainforest zone lies between latitudes 5.91˚and 9.29˚N and longitude 2.79˚and 6.11˚E [29]. Temperature ranges between 21˚C and 34˚C while the annual rainfall ranges between 1500 mm and 3000 mm. The vegetation consists of fresh water swamp and mangrove forest at the belt, the lowland forest, and secondary forest.

Sampling procedure
A total of 569 adult (8 months old) Nigerian indigenous Guinea fowls comprising 109 birds (27 males and 82 females) from Southern Guinea Savanna zone, 270 birds (80 males and 190 females) from Sudano-Sahelian zone and 190 birds (51 males and 139 females)from Tropical Rainforest zone were used in the study. The indigenous birds were randomly sampled in smallholder rural farmers flocks and managed under the traditional low-input settings.
Multistage sampling procedure was purposively and randomly adopted in the selection of States, Local Government Areas (LGAs), villages and Guinea fowl keepers in each agro-ecological zone. States, LGAs, and villages were purposively selected based on the knowledge of the availability of Guinea fowls in the communities as provided by the local Extension Agents and Community Heads. The number of sampling locations varied with 4 LGAs and 11 villages in the Southern Guinea Savanna, 5 LGAs and 15 villages in the Sudano-Sahelian, and 4 LGAs and 13 villages in the Tropical Rainforest. Based on willingness to participate in the research, eleven individuals were then randomly selected from each village making a total of 429 households (n = 121, 165 and 143 for Southern Guinea Savanna, Sudano-Sahelian and Tropical Rainforest, respectively).

Data collection
Data collection was done in the rainy season (April to June, 2020). Morphologically distinct Guinea fowls were identified using phenotypic traits based on the standard descriptors by FAO [14], AU-IBAR [30] and the colour chart of Guinea fowl by GFIA [31]. The sexes were distinguished through visualisation of the vent and the use of helmet shape as well as wattle size and

PLOS ONE
shape [25]. Eleven qualitative physical parameters such as plumage colour, skin colour, shank colour, eye colour, earlobe colour, helmet colour, helmet shape, wattle possession, wattle size, wattle shape and skeletal structure were used to characterize the Guinea fowls morphologically. For quantitative (biometric) description, the following body parts were measured: Body weight (kg): The live weight of the Guinea fowl; Head length (cm): Taken between the most protruding point of the occipital and the frontal (lacrimal) bone; Head thickness (cm): Head thickness measured as the circumference at the middle of the head; Helmet length (cm): Measured as the distance between the base of the head to the tip of the helmet; Helmet width (cm): Measured as the distance between the broadest part of the helmet; Wattle length (cm): Taken as the distance between the base of the beak and the tip of the wattle; Wattle width (cm): Measured as the distance between the broadest part of the wattle; Neck length (cm): Distance between the occipital condyle and the cephalic borders of the coracoids; Neck circumference (cm): Taken at the widest point of the neck; Wing length (cm): Taken from the shoulder joint to the extremity of the terminal phalanx, digit III; Wing span (cm): Distance between the two wings when stretched out; Body length (cm): The distance from the first cervical vertebra (atlas) to the posterior end of the ischium; Trunk length (cm): The distance between shoulder joint and posterior edge of the ischium; Keel length (cm): Keel length (sternum or breast bone) measured from the anterior point of the keel to the posterior end; Chest circumference (cm): Taken as the circumference of the body around the breast region; Thigh length (cm): Distance between the hock joint and the pelvic joint; Shank length (cm): Measured as the distance between the foot pad and the hock joint; and Shank thickness (cm): Measured as the circumference at the middle or widest part of the shank.
Also, the following morphological indices were estimated [32]: Massiveness: The ratio of live body weight to trunk length x 100; Compactness: The ratio of chest circumference to trunk length x 100; Long-leggedness: The ratio of shank length to body length x 100; and Condition index: The ratio of live body weight to wing length × 100.
The weight measurement was taken using a hanging digital scale (WeiHeng Brand), the width measurements were taken using a vernier caliper (0.01 mm precision) while the length and circumference measurements were taken using a flexible tape measure.

Statistical analysis
Descriptive statistics. Descriptive statistics were computed to determine the frequencies of the qualitative physical traits. Where statistical significant differences in the frequencies were obtained at agro-ecological and sex levels, they were assessed using the non-parametric Kruskal-Wallis H test followed by the Mann-Whitney U test for post hoc separation [33] of IBM-SPSS software [34]. This approach was adopted as a result of the small and unequal sample sizes of phenotypic groups including non-normality of the data distribution.
Correspondence analysis. Multiple correspondence analysis (MCA) was used to establish the relationships between the qualitative physical traits using JMP 16 [35] statistical software. In order to run the MCA, the input data (qualitative physical traits including their classes) were saved in IBM-SPSS software and opened under the JMP input file platform. Then, MCA was selected under multivariate methods. Preliminary analysis revealed that wattle possession and skeleton structure had zero variance and were excluded from the MCA.
Univariate analysis. Biometric traits and morphological indices were tested for normality with the Shapiro-Wilk's test (P<0.05) and by visual inspection of the histograms. Levene's test was used to confirm homogeneity of variances (P>0.05) as decribed by Brown et al. [36]. Due to small and unequal sizes, low male-female ratio and non-normality of the distribution of the data, the non-parametric Kruskal-Wallis H test was performed to compare mean ranks of biometric traits and morphological indices based on agro-ecology, sex, and sexes within each agro-ecology. In the case of significant Kruskal-Wallis H test, Dunn-Bonferroni test (agroecology) and Mann-Whitney U test (sex, and sexes within each agro-ecology) were used for pairwise comparisons of mean ranks.
Stepwise canonical discriminant analysis. Canonical discriminant analysis [37] option of IBM-SPSS [34] statistical software was applied to classify birds in the three agro-ecological zones based on quantitative traits. In the analysis, all the eighteen biometric traits and four morphological indices (covariates) were entered in a stepwise fashion as explanatory variables to establish and outline population clusters [38] based on agro-ecology. F-to-remove statistics was the criterion for variables' selection while multicollinearity was detected among the variables in the discriminant function using tolerance statistics. The ability of this discriminant model to identify birds in the Southern Guinea Savanna, Sudano-Sahelian and Tropical Rainforest zones was indicated as the percentage of individuals correctly classified from the sample that generated the model. The accuracy of the classification was evaluated using split-sample validation (cross-validation).
Categorical principal component analysis. Categorical principal component analysis (CATPCA) procedure was employed to explore hidden relationships among the qualitative physical traits (with the exception of wattle possession and skeleton structure due to zero variance), biometric traits and morphological as described by Martin-Collado et al. [39]. This was to allow for appropriate grouping of the guinea fowls based on agro-ecology and sex. The PCs were extracted based on Eigenvalues greater than 1 criterion. The convergence was 0.00001 with maximum iterations of 100. The PC matrix was rotated using the varimax criterion with Kaiser Normalization to facilitate easy interpretation of the analysis. The reliability of the PCA was tested using Chronbach's alpha using IBM-SPSS [34].
Decision trees. Chi-square automatic interaction detection (CHAID) and Exhaustive CHAID algorithms were employed to assign the birds into agro-ecological zones using the qualitative physical traits (with the exception of wattle possession and skeleton structure due to zero variance), biometric traits and morphological indices as the predictor variables. CHAID is a tree-based model with merging, partitioning and stopping stages that recursively uses multi-way splitting procedures to form homogenous subsets using Bonferroni adjustment until the least differences between the predicted and actual values in a response variable are obtained [40]. It produces terminal nodes and finds the best possible variable or factor to split the node into two child nodes. The Exhaustive CHAID, as a modification of CHAID algorithm, applies a more detailed merging and testing of predictor variables [41]. The accuracy of CHAID and Exhaustive CHAID models was obtained from the percentage of individuals correctly classified in each agro-ecological zone. The predictive performance of each model was assessed using the goodness-of-fit criteria [40]. The most predictive model estimates the highest values in correlation coefficient (r), coefficient of determination (R 2 ) and Adj R 2 , and the lowest values in relative approximation error (RAE), mean absolute error (MAE), standard deviation ratio (SDratio), root mean square error (RMSE) and the coefficient of variation (CV, %), respectively. IBM-SPSS [34] software was also used for the Decision Trees' analysis

Distribution of the qualitative traits
The frequency distribution of the colour traits of indigenous helmeted Guinea fowl is shown in Table 1. Agro-ecology significantly affected (P<0.05; P<0.01) all the six traits investigated. No definite pattern of variation in each class of the colour traits was observed among the three agro-ecological zones. Generally, the most frequent colour phenotype of helmeted Guinea fowl in Nigeria had pearl plumage colour (54.0%), pale red skin colour (94.2%), black shank colour (68.7%), brown eye colour (49.7%), white earlobe colour (54.8%) and helmet colour (72.6%). However, sex did not influence (P>0.05) all the six colour traits.
The frequencies of helmet shape and wattle size were significantly affected by agro-ecology (P<0.01) ( Table 2). While most of the birds had single helmet shape (50.8%), which appeared to be more in the Sudano-Sahelian and Tropical Rainforest zones, wattle size did not follow a definite pattern. All the birds in the three agro-ecologies had wattle and were skeletally normal (P>0.01). However, sex had a significant effect (P<0.01) on helmet shape (where more females were single), wattle size (where that of males appeared larger), and wattle shape (where more females carried theirs flat).

Biplot of the multiple correspondence analysis
The MCA revealed the association between the qualitative physical traits and agro-ecological zones in two dimensions (Fig 2). The first dimension was high and represented 93.2% of the deviation from independence while the second dimension signified 6.8% of the total variation based on the inertia. The agro-ecological zones were not clustered perfectly (as revealed bythe low inertia values of 0.168 and 0.012) considering the intermingling of some qualitative physical traits. This was more noticeable between birds in the Sudano-Sahelian and Tropical Rainforest zones. Therefore, discrimination of the traits appears very weak. However, on the right hand side of the biplot, peach black, orange and pale pink shank colour, dark skin colour, and red and slanted backward helmet seemed to be more associated with the Southern Guinea Savanna zone.

The fixed effect of agro-ecology on biometric traits and morphological indices
The results of the univariate analysis revealed significant effect (P<0.05) of agro-ecology on the biometric traits and morphological indices of the guinea fowls [Medians (means in parentheses)] (Table 3). Overall, birds from the Southern Guinea Savanna zone had significantly higher values (P<0.05) for most zoometrical traits compared to their Sudano-Sahelian and Tropical Rainforest counterparts. However, the former and the latter were similar (P>0.05) in Table 1

Medians (means in parentheses) of biometric traits and morphological indices of indigenous helmeted Guinea fowls based on agro-ecology.
Agro-ecology

The fixed effect of sexes within agro-ecologies on biometric traits and morphological indices
The effect of sexes within agro-ecologies had significant effect (P<0.05) on some biometric traits and morphological indices in two out of the three agro-ecological zones [Medians (means in parentheses)] (

PLOS ONE
long-leggedness, neck length, thigh length and wattle length were the significant (P<0.001) parameters of importance to separate birds in the Southern Guinea Savanna, Sudano-Sahelian and Tropical Rainforest zones. However, there was considerable spatial intermixing of the biometric traits largely observed between birds in the Sudano-Sahelian and Tropical Rainforest zones (Fig 3). The predicted group membership of the three agro-ecological zones is shown in Table 7. The classification results showed that 88.1, 51.9 and 55.8% of birds in the Southern Guinea Savanna, Sudano-Sahelian and Tropical Rainforest zones, respectively were correctly assigned to their distinct groups. The three respective group cases were 57.1% cross-validated.

Contributions to variation and loadings of variables on the principal components
The result of CATPCA revealed the extraction of two principal components (PCs) which explained 42.1% of the variation in the dataset ( Table 8

Decision trees of the data mining
The tree diagram of the CHAID algorithm is shown in Fig 5. Seven terminal nodes (Nodes 1, 2, 3, 5, 6, 7 and 8) were formed. The root node (Node 0) showed the descriptive statistics of the birds in the three agro-ecological zones. The Chi-squared-based branch and node distribution   (Fig 6). Here, wing length (>18.10 cm) was also the best single discriminant variable (P<0.001) to distinguish birds in the three agro-ecological zones. In contrast to what was obtained under CHAID, body length and eye colour were the two additional variables to differentiate the populations. Wing length (14.80-15.50 cm), body length (< = 35.00 cm) and eye colour permitted a better separation of the Sudano-Sahelian from Tropical Rainforest birds. Unlike what was observed in CHAID, birds from the former had mostly brown, white and pink eye colours while the later were characterized by black as well as bluish eye colour. In this model, 52.3, 77.4 and 37.9% of birds in the Southern Guinea Savanna, Sudano-Sahelian and Tropical Rainforest zones, respectively, were correctly assigned to their distinct agro-ecological zone with an average accuracy rate of 59.4% (Table 10). The r, R 2 , Adj R 2 , RAE, MAE, SDratio, RMSE and CV (%) values were 0.520, 0.270, 0.268, 0.282, 0.009, 0.896, 0.637, and 29.765, respectively.

Discussion
Phenotypic variation of local animal resources indicates a genetic diversity that may be worth conserving for future uses while better understanding of the external features helps to facilitate

PLOS ONE
the implementation of conservation policies aimed to ensure local resources survival [15]. Morphometric and phaneroptic approaches may be fundamental in the management of poultry, considering the fact that they are fast and economically profitable [37]. The preponderance of more female birds in the present study could be attributed to the fact that smallholder poultry farmers normally keep more hens for the purpose of procreation, whereas the cocks are mostly slaughtered for consumption or sold to generate family income. We observed four major plumage colours (Pearl, Black, Brown and Pied). The varying colour patterns could be an indication that there are no pure genotypes of Guinea fowl in Nigeria as there are no records of selective breeding of the indigenous stock birds. However, the colour patterns were somehow different from the dominant Pearl, Lavender, Black and White variations earlier reported in the country [42,43]. The slight variation may be occasioned by sampling coverage.
In a similar study in Ghana, Agbolosu et al. [21] found that the predominant plumage colour was pearl grey colour (43.7%), whereas Traore et al. [22] reported pied plumage colour (42.76%) as the most frequent colour in the provinces of Burkina Faso. The Nigerian birds shared brown eye colour (57.0%) with those of Atakora (Mountainous) dry savannah zone in Togo [44] and black shank colour with those of Kenya (95.6%) [16], Sudanian and Sudano-Guinean zones in Benin [45]. Colour polymorphism defies evolutionary expectations as a

PLOS ONE
single species may maintain a striking phenotypic variation [46]. The present variant phenotypes may be due to polymorphism [47] and might have evolved in local Guinea fowls as adaptive measures for survival under varied environmental conditions. According to Getachew et al. [48], sustainable livestock production in the tropics requires adaptive genotypes which can withstand the undesirable effects of climate change and ensure optimal performance of the birds. In another study on a different species, Nigenda-Morales et al. [49] reported that the overall fitness of individuals in their environments may be affected by colour while Gong et al. [50] considered colour variation as an environmental indicator, which provides clues for the study of population genetics and biogeography. The preponderance of Pearl plumage colour in our study may also be attributed to farmers' preference, which is congruous with the submission of Vignal et al. [4] that prevalence of a particular colour could be attached to socialcultural value without any proven relationship with a biological function. This was buttressed by the report of González Ariza et al. [37] that certain qualitative physical traits may be associated with consumers' trends and their cultural preferences. Our findings on helmet shape are

PLOS ONE
in agreement with the report on indigenous Guinea fowls in Ghanian where single shape (42.70%) predominated. The current observation on helmet shape where more females exhibited single shape is congruous withthe submission of Angst et al. [51] that females have bony helmet more compact dorsoventrally while the males have taller helmet, with a more complex shape including curvature of the posterior part along the dorsoventral axis. Similarly, Agbolosu et al. [21] reported that helmet shape is more pronounced in males than females. The observation on wattle is in consonance with the findings of Umosen et al. [52] who stated that, on the average, females had small wattle which was mostly flat. In order to ascertain the genetic purity of the birds, the MCA result did not give a perfect clustering of the birds as phenotypic homogeneity of the Guinea fowl populations was evident in Sudano-Sahelian and Tropical Rainforest birds. This is in spite of the wide geographical distance and varying environmental conditions between the two zones. This suggests that colour traits alone might not be enough to distinguish between the three agro-ecological zones. Similar submission was made by Traore et al. [22] where, in spite of the enormous environmental differences, there was morphological homogeneity in qualitative traits in Guinea fowls in Burkina Faso. Brown et al. [36] also observed limited phenotypic and genetic diversity in local Guinea fowls in northern Ghana.
Univariate analysis revealed significant differences among zones for most biometric traits and calculated body indices, suggesting the possible influence of these zones on the evolutionary adaptation of the Guinea fowl population. However, there was no clear cut pattern in the biometric traits and morphological indices especially of the Sudano-Sahelian and Tropical Rainforest birds. The body weights of the present study are comparable to the 1.40 kg reported by Orounladji et al. [45] for indigenous Guinea fowls in a Sudanian zone in Benin. They are however, higher than the range 1.08-1.33 kg reported for adult Guinea fowl (Numida meleagris) in a humid zone of southern Nigeria [53] and 1.275 kg obtained in Zimbabwe [54]. Nevertheless, the indigenous birds are smaller in size when compared to their exotic counterparts. While Agwunobi and Ekpenyong [55] obtained a live weight of 1.5 kg for 'Golden Sovereign' Guinea fowl broiler strain under tropical conditions of Nigeria, Batkowska et al. [56] found a range of 2166 ± 42.5-2291 ± 46.9 kg for French commercial set. The differences may be attributed to genetics, age, physiological stage, location and management systems employed by the keepers. According to Ahiagbe et al. [57], genetic make-up and management practices could affect the growth traits of Guinea fowls. Exotic Guinea fowls are products of many years of robust selection and breeding [58,59]. Therefore, it is possible that crossbreeding between the indigenous and exotic will result in birds of high genetic superiority in terms of meat yield and quality, egg production and adaptation. Sexual dimorphism provides insight into the sexualand natural-selection pressures being experienced by male and female animals of different species [60]. At inter-population level, especially with some biometric traits, sexual dimorphism in the present study favoured male animals. This concurs with the established literature that males generally possess larger body sizes than females in normal sexual size dimorphism in birds [61]. The differential rate and duration of growth by the sexes may be responsible for the PLOS ONE present observations. Also, high rate of breeding in the populations could be another contributing factor to sexually dimorphic traits [62], as the birds have not been selected for the purpose of classical breeding. As obtained in the current study, Dudusola et al. [53] found male dominance in thigh length, body length, wing length, wing span, wattle length and chest circumference in Nigeria while Brown et al. [36] reported longer body and shank length including wingspan in indigenous Guinea fowl in Ghana. In a related study on domestic chicken, Toalombo Vargas et al. [63] reported longer thigh length in male birds. The canonical discriminant analysis showed high level of admixture especially between the Sudano-Sahelian and Tropical Rainforest populations. It could, therefore, be reported that the Guinea fowls in Nigeria are unselected and largely of mixed populations. Northern Nigeria is the traditional home of indigenous helmeted guinea fowls in the country [64]. Considering the geographical proximity of the Southern Guinea Savanna and Sudano-Sahelian zones, one would have expected considerable intermixing of the guinea fowl populations. However, the reverse was observed in the present study as the intermingling between the birds in the Sudano-Sahelian and Tropical Rainforest zones was higher which could partly be due to transhumance especially by herders. The herders (mainly cattle rearers) from the northern parts of the country do move to the southern parts in search of natural pastures during the dry season. When they do so, they tend to carry along all their animals to their new locations. In that process, there is the possibility of exchange of birds between the settlers and their hosts. Such livestock mobility, which is seen as a means to an end [65] could have shaped poultry distribution pattern. Suffice to say that the guinea fowl (Numida meleagris) population of Tropical Rainforest is an ecotype of the Sudano-Sahelian; which is quite different from Numida ptilorhycha that is indigenous to the deciduous rain forest zone of southern Nigeria [66]. This assertion is corroborated by the reports of Ayorinde [67] and Obike et al. [68] who observed that Numida meleagris, domiciled in the north was spreading to other smallholder farming areas. In a related study, Whannou et al. [69] submitted that the mobility of herders could engender genetic introgression, thereby affecting animal genetic diversity. Another possible factor that could have contributed to the genetic erosion is inter-regional trade. It appears such live animal trade seemed to be more between livestock marketers in the Tropical Rainforest and Sudano-Sahelian zones than their Southern Guinea Savanna counterparts. According to Benton et al. [70], market dynamics in one location could drive biodiversity-damaging practices in other locations. In another study, Valerio et al. [71] highlighted the relevance of cross-border ties suggesting that markets play distinct structural roles in understanding animal movement patterns.
The results of CATPCA showed that some levels of separation of the Guinea fowls can be obtained based on agro-ecology which was more associated with the second principal component. The body parameters of importance in this component are wing length, skin colour, neck length, head length, wing span, eye colour, shank thickness, helmet colour, helmet shape, earlobe colour, wattle size, plumage colour, keel length, shank colour and trunk length. These parameters describe more of shape and colour of the guinea fowls. However, these differences in biometric traits and morphological indices based on agro-ecology were weak due to the fact that the second principal component could only account for 15.0% of the total variation. The use of CATPCA in assigning birds to their genetic groups had earlier been reported [33].
The decision tree results revealed that the guinea fowls from the Southern Guinea Savanna, Sudano-Sahelian and Tropical Rainforest zones could to a relative extent be separated using wing length, body length, and eye colour. However, the average accuracy rate of 58.0% (CHAID) and 59.4% (Exhaustive CHAID) obtained in this study indicated that 42 and 31.6% of the birds were wrongly classified. The implication of this is that there is a form of intermixing of the birds in the three agro-ecological zones. Both wing and body lengths are skeletal parameters that are not influenced by body condition, thereby providing good estimates of overall body size of the birds. It is possible that both traits are under similar selection pressure [72]. The importance of morphometric traits in population stratification has also been stressed in other avian species [73,74].
When all the algorithms used in this study are jointly considered, it could be said that the Guinea fowls from the Southern Guinea Savanna, Sudano-Sahelian and Tropical Rainforest zones of Nigeria are more homogeneous than heterogeneous in terms of the investigated qualitative physical traits, biometric traits and morphological indices. The biological implication of this is that elite birds from the three agro-ecological zones could be selected for the purpose of pure breeding or crossbreeding with their more productive exotic counterparts. This is beneficial considering the fact that the existence of several varieties of Guinea fowls on farms does not encourage their genetic conservation and improvement [75]. Our present findings are similar to the report of Traore et al. [76], where Guinea fowls in Burkina Faso were highly intermingled, suggesting that differences in biometric and qualitative physical traits were not related to geography. In a related study, Etienne et al. [77] reported that local Guinea fowls in Côte d'Ivoire exhibited less phenotypic diversity. In another study, it was found that Guinea fowls in northern Togo belonged to a single indigenous population [78].

Conclusion
The qualitative physical traits of Nigerian Guinea fowls predominantly were affected by agroecology. However, there was no clear cut variation and distribution pattern across the three agro-ecological zones. Although the indigenous birds generally were of low body weights, those in the Southern Guinea Savanna zone were more compact while their counterparts in the Sudano-Sahelian and Tropical Rainforest zones had longer legs than body, and better condition index. Small body size could be part of the animals' adaptation for survival under the low-inputs tropical environment. The superiority of male birds to their female counterparts could be attributed to sexual dimorphism. The clustering pattern of the traits based on MCA and canonical discriminant analysis especially between the Sudano-Sahelian and Tropical Rainforest birds revealed high level of admixture, although the bird populations to an extent could be distinguished using wing length, body length, and eye colour. Overall, it could be said that the guinea fowls from the three agro-ecological zones exhibited less phenotypic diversity, and belonged to a single indigenous population. However, there is a need for further genomic studies to consolidate the present findings, and pave the way for policy decisions geared towards effective management, conservation and genetic improvement of the indigenous birds. The anticipated benefits include the development of hybrid improved Guinea fowls for the empowerment of women and youth including improvement in food security and livelihoods.